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^vj ■ Abstract. We review the main properties of shell models for magnetohydrody- 

namic (MHD) turbulence. After a brief account on shell models with nearest neigh- 
'~{ , bour interactions, the paper focuses on the most recent results concerning dynam- 

ical properties and intermittency of a model which is a generalization to MHD 

^^ ' of the Gledzer-Yamada-Okhitani (GOY) model for hydrodynamic. Applications to 

astrophysical problems are also discussed. 

Q . 1 Introduction 

u 

f-i . Shell models are dynamical systems (ordinary differential equations) rep- 

resenting a simplified version of the spectral Navier-Stokes or MHD equa- 
tions for turbulence. They were originally introduced and developed bydfg 
Obukhov |]1|, Desnyansky and Novikov [|| and Gledzcr ||] in hydrodynamic 
turbulence and constitute nowdays a consistent and relevant alternative ap- 
^ ' proach to the analytical and numerical study of fully developed turbulence 

CN . (see [H for a complete review). 

f^ ' Shell models are built up by dividing wave-vector space (k-space) in a 

discrete number of shells whose radii grow exponentially like k„ = fcgA", 

V— s ' (A > 1), n = 1, 2, . . . , iV. Each shell is assigned a scalar dynamic variable, 

Unit), (real or complex) which takes into account the averaged effects of 
velocity modes between fc„ and fc„+i. The equation for M„(i) is then written 
in the formex 

^ = KCn + Dn + Fn (1) 

at 

where knCn, Dn and F^ are respectively quadratic nonlinear coupling terms 

(involving nearest or next-nearest shell interactions), dissipation terms and 

forcing terms, the last generally restricted to the first shells. Nonlinear terms 

H , are chosen to satisfy scale- invariance and conservation of ideal invariants. The 

main advantage shell models offer is that they can be investigated by means 
of rather easy numerical simulations at very high Reynolds (Re) numbers. 
The degrees of freedom of a shell model are N ^ In Re, to be compared with 
N ^ Re^/'^ for a three dimensional hydrodynamic turbulence following the 
Kolmogorov scaling. 

The paper is organized as follows. In section 2 shell models with nearest 
neighbour interactions are briefly reviewed. In section 3 equations for MHD 
models with nearest and next-nearest neighbour interactions are presented 
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and conservations laws for the ideal case are discussed. Section 4 is devoted to 
dynamo action in shell models and section 5 to spectral properties in forced 
stationary state and intermittency In section 6 conclusions are drawn and a 
brief mention to astrophysical applications is made. 



2 Models with nearest neighbour interactions 

The simplest hydrodynamic shell model is the Obukhov-Novikov model, 
which is a linear superposition of the Obukhov equation |l| and the Novikov 
equation M. The model involves real variables u„(t) and conserves the en- 
ergy 1/2 ^„^]^ M^ in absence of forcing and dissipation. It does not conserve 
phase space volume nor other quadratic invariants exist. The extension of 
the Obukhov-Novikov model to MHD is due to Gloaguen et al. [p| . We write 
down the equations for clarity (u„ and 5„ represent respectively the velocity 
and the magnetic field in dimensionless units) 

— ^ = -vkl_Un + a {knul_-^ - kn+lUnUn+1 - fcn&^_i + fcn+l^«6n+l) (2) 



dt 

db„ 
dt 



+13 {knUn-lUn - kn+iU^^^ - fc„6„6„_i + kn+ib^^j) 
= -V^lbn + akn+1 {Un+lbn - U„6„+l) + Pk^iu^bn-l - Un-lbn) (3) 



Here ly is the kinematic viscosity, rj is the magnetic diffusivity, a and P are two 
arbitrary coupling coefficients. The ideal invariants of the system are the total 
energy, 1/2 J2n=i ( "n + ^n ) ^^'^ the cross-correlation, J2n=i '^nbn which are 
two ideal invariants of the MHD equations 1^ . When written in terms of the 
Elsasser variables Z^ = u„ + bn-, Z~ = m„ — 6„, the equations assume a sim- 
metric form and the conservation of the two previous invariants is equivalently 
expressed as the conservation of the pseudo-energies E^ = (1/4) X]n=i ^n ^• 
It is remarkable to note that, unlike the hydrodynamic model, the MHD ver- 
sion satisfies a Liouville theorem J2n=i 9{dZ^/dt)/dZ^ — 0, imphng phase- 
space volume conservation. The MHD equations conserve a third ideal invari- 
ant which is the magnetic helicity in three dimensions (3D) and the mean 
square potential in two dimensions (2D), but no further ideal quadratic in- 
variant can be imposed to this shell model. 

A detailed bifurcation analysis for a three-mode system was performed in 
[|j for different values of a and /3. The low Reynolds (kinetic and magnetic) 
numbers, used as control parameters, allowed to identify a great variety of re- 
gions in the parameter space. Turbulence was investigated with a nine-mode 
system which produces an inertial range with spectra following approximately 
the Kolmogorov scaling E{k) ^ k~^^^. Temporal intermittency was also ob- 
served and then reconsidered in more details by Carbonc m who calculated 
the scaling exponents of the structure functions for the Elsasser variables 
and for the pseudo-energy transfer rates, showing consistency with the usual 
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multifractal theory. Other interesting MHD phenomena were also observed in 
[pi such as dynamo effect and the growth of correlation between velocity and 
magnetic field in an unforced simulation. These phenomena will be treated 
in more details in the next paragraphs. 

The complex version of (0) and (H) was thoroughly investigated by Biskamp 
[B| . The complex model allows to include the Alfven effect M , fTq] , [nil , that 
is the interaction of a constant large scale magnetic field with small scale 
turbulent eddies. The main consequence of this effect should be a reduction 
of the spectral energy transfer rate and a consequent change of the spectra 
from the Kolmogorov scaling, E{k) ^ k~^'^, to the Iroshnikov-Kraichnan 
one, E{k) ^ fc"'^/^. In this paper the Alfven effect will not be furtherly 
treated. The reader is referred to [|8| for a complete discussion concerning the 
inclusion of Alfvenic terms in shell models. 

3 Models with nearest and next-nearest interactions 

Shell models with nearest and next-nearest neighbour interactions were in- 
troduced by Gledzer H. In particular the so called GOY (Gledzer-Yamada- 
Ohkitani) model has been extensively both numerically and analitically inves- 
tigated |13| , |13| , [njl . The GOY model allows to conserve another quadratic 
invariant besides energy which was identified with the kinetic helicity u5\. 
A generalization of the GOY model to MHD can be found in Biskamp^ |. 
All the parameters of the model are now fixed by imposing the conservation 
of another quadratic invariant that can be chosen to distinguish between a 
3D and a 2D model. A more refined version was then considered by Frick 
and Sokoloff ^^ to take into account the fact that the magnetic helicity is 
a quantity not positive definite. The situation can be summarized as follows 



Let us consider the following set of equations (m„ and 6„ are now complex 
variables representing the velocity and the magnetic field in dimensionless 
units) 

—^ = -lyk^Un - v'k^^'^Un + ikn\{Un+lUn+2 - bn+lbn+2) 

5 1 — 5 \* 

--(w„_lW„+l - fo„-lfcn+l) T5-(m«-2M«-1 - ^«-2fon-l) j + fn (4) 

—n- = ^Vklbn + iknUl - 5 - S,n){Un+lbn+2 - &n+lM„+2) 

6 1 — S 'I * 

+ -Yi'^n~lbn+l -bn-lUn+l) -\ TJ^ {Un-2bn-l - bn-2Un-l) > + 9,i (5) 

or, in terms of the complex Elsasser variables Z^{t) — Vn{t) ± bn{t), partic- 
ularly useful in some solar-wind applications, 

dZ^ v' v' 

" =-v+klzt-v-klzi--k-^z^,--ki^z-^ikrj:t + !n (6) 



^ ,„^.^^ ^ .^^^^ ,^^ ^^^ 
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where 



2 '^n+l^n+2 + 2 ^ri+l^n+2 



2A "+1 n-1 2A n+1 n- 

(5„i - (5 „ 1 _^ 2 - (5 - (5„ 



2_\2 -^n-l^n-2 2A2 ^n-l^n-2 } C^) 

Here v = [v ± ?7)/2, being i' the kinematic viscosity and 77 the resistivity, 
—v'k~'^Un, eq. (Q), is a drag term specific to 2D cases (see below), f^ — 
(fn ign)/'^ are external driving forces, S and 6m are real coupling coefficients 
to be determined. In the inviscid unforced limit, equations (ra) conserve both 
pseudoenergies E^{t) = (1/4) X^n l^rt (^)l for any value of S and Sm (the 
sum is extended to all the shells), which corresponds to the conservation of 
both the total energy E = E+ + E- = (1/2) ^„(|w„(i)|^ + \bn{t)\^) and 
the cross-helicity he = E^ — E~ = J^n^^i'^^bn)- As far as the third ideal 
invariant is concerned, we can define a generalized quantity as 



N 

Hi-\t) = Y^{sign{5-l)r^-^^^ (8) 



r(«)^+^ - "^ /..;„„ /^ T\\n\bn(f)? 



ha 

whose conservation implies 5=1 — A^", 5„i = A^"/(l + A~") for 6 < 1, 
< ,5„ < 1 and ^ = 1 + A-", (5,„ = -A-"/(l - A"") for 6 > 1, d,n < 0, 
6m > 1- Thus two classes of MHD GOY models can be defined with respect 
to the values of S: 3D-like models for ^ < 1, where Hg is not positive def- 
inite and represents a generalized magnetic helicity; 2D-like models where 
6 > I and H^' is positive definite. This situation strongly resembles what 
happens in the hydrodynamic case where 2D-like {5 > 1) and 3D-like {6 < 1) 
models are conventionally distinguished with respect to a second generalized 
conserved quantity Hj^'{t) = J2n (sign((5 — 1))"^;" |t'ra(i)| ■ Here the 3D and 
2D cases are recovered for a = 1, 2 where the ideal invariants are identified 
respectively with kinetic helicity and enstrophy. It should be noted that, al- 
though the hydrodynamic invariants are not conserved in the magnetic case, 
the equations which link a and S are exactly the same for hydrodynamic and 
MHD models. Thus, once fixed a and i5, it is a simple matter to find out which 
GOY model the MHD GOY one reduces to when bn = [|7|. To summarize 
we have that (with A = 2) the model introduced in jl^ for the 3D case will be 
called, hence on, 3D MHD GOY model or simply 3D model. It is recovered 
for a = 1, 5 = 1/2, Sm = 1/3 and reduces to the usual 3D GOY model for 
bn = 0. The Biskamp's 3D model H is actually a 2D-like model and will be 
called pseudo 3D model. It is obtained for a = 1, S = 3/2, 6m — —^ and 
reduces to a 2D-like GOY model that conserves a quantity which has the 
same dimensions as kinetic helicity but is positive definite. The 2D models 
introduced in S and in iQ coincide, they are recovered for a = 2, 5 = 5/4, 
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6m = — 1/3 and reduce to the usual 2D GOY model for 6„ = 0. In the fol- 
lowing the properties of the 3D model will be mainly investigated. 



4 Dynamo action in MHD shell models 

The problem of magnetic dynamo, that is the amplification of a seed of 
magnetic field and its maintenance against the losses of dissipation in an 
electrically conducting flow, is of great interest by itself and for astrophys- 
ical applications (see for example |Q for an excellent introduction to the 
problem). Shell models offer the opportunity to test with relative simplicity 
whether a small value of the magnetic field can grow in absence of forcing 
terms on the magnetic field. Previous considerations about dynamo action 
in shell models can be found in g . In that case numerical study of bifurca- 
tions in the three-mode system revealed instabilities of kinetic fixed points to 
magnetic ones or magnetic chaos. The existence of a sort of dynamo effect in 
MHD GOY models was put forward by Frick and Sokoloff iQ . The authors 
investigate the problem of the magnetic field generation in a free-decaying 
turbulence, thus showing that: 1) in the 3-D case magnetic energy grows and 
reaches a value comparable with the kinetic one, in a way that the magnetic 
field growth is unbounded in the kinematic case; 2) in the 2D case magnetic 
energy slowly decays in the nonlinear as well as in the kinematic case. These 
results have been interpreted as a 3D "turbulent dynamo effect" and seem 
to be in agreement with well-known results by which dynamo effect is not 
possible in two dimensions p9| . The problem was then reexamined in p^ in 
a forced situation looking at a comparison between the 3D MHD GOY model 
and the pseudo 3D model. 

Starting from a well developed turbulent velocity field, a seed of magnetic 
field is injected and the growth of the magnetic spectra monitored. System 
is forced on the shell n = 4 (fco = 1), setting f^ = f^ = {l + i) 10~^, which 
corresponds to only inject kinetic energy at large scales. Method of integration 
is a modified fourth order Runge-Kutta scheme. In fig. |l| we plot logio( l^"l^ ) 
and log;^o( kn|^ ) versus log^g fc„ for the 3D model. Angular brackets ( ) stand 
for time averages. It can be seen that the magnetic energy grows rapidly 
in time and forms a spectrum where the amplitude of the various modes 
is, at small scale, of the same order as the kinetic energy spectrum. (The 
subsequent evolution of magnetic and kinetic spectra will be considered in the 
next section). The spectral index is close to k~'^/'^ which is compatible with a 
Kolmogorov scaling of the second order structure function. For a comparison 
we integrated the pseudo 3D model and it can be seen (fig. ||) that a magnetic 
spectrum is formed, but it slowly decays in time. Notice that, because of the 
smallness of fe„, its back-reaction on the velocity field is negligible, thus the 
kinematic part of the model evolves independently from the magnetic one. 
Now the scaling |w,ip ~ fc„ follows, as a cascade of generalized enstrophy 




Fig. 1. 3D model: logj^Q( |unp ) (diamonds) and logj^o( l^nf ) (lines) versus logj^Q fc„. 
The averages of |6„P are made over intervals of 3 large scale turnover times. Time 
proceeds upwards. The kinetic spectrum is averaged over 30 large scale turnover 



times. The straight line has slope 



-2/3. Parameters used: N=24, u = t] — 10 



is expected for 2D-like hydrodynamic GOY models when a < 2 (see |20 for 
details). The question now arises whether it is correct the interpretation of the 
growth of the magnetic field in the 3D model as the corresponding dynamo 
effect expected in the real 3D magnetohydrodynamics. First of all it should 
be noted that in the kinematic case an analogy with the vorticity equation 
predicts the following relations between velocity and magnetic energy spectra 
[|j: Ivn]"^ ~ fc"", |6„1^ ~ fc^~°, so that if a = 2/3 it follows a magnetic energy 
spectrum growing with k. The kinematic case corresponds to the first stage of 
growth of our simulation where this behaviour is sometimes visible, at least 
qualitatively. Note however that the averages are made on very small time 
intervals because of the rapid growth of the magnetic energy. A similar, much 
more pronounced behaviour is found for the pseudo 3D model as well. 

Let us stress that the sign of the third ideal invariant seems to play a 
crucial role as far as the growth of small magnetic fields is concerned. In ef- 
fect this sort of dynamo effect can also be considered under a different point 
of view. Let us consider the ideal evolution of the model dZi^/dt = iknT^* ■ 
We can build up the phase space S of dimension D = AN, by using the 
Elsasser variables as axes, so that a point in S represents the system at a 
given time. A careful analysis of (p) shows that there exist some subspaces 
I C S oi dimension D = 2N which remain invariant under the time evo- 
lution pl|]. More formally, let y{0) — (ti„,5„) be a set of initial conditions 
such that y(0) € I, I is time invariant if the flow T*, representing the time 
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Fig. 2. Pseudo 3D model: logio( \vn\^ ) [diamonds) and logj^Q( \bn\^ ) (lines) versus 
logj^o ^n- Averages are made over intervals of 100 large scale turnover times. Time 
proceeds downwards. The kinetic spectrum is only shown for the last interval. The 
straight line has slope —4/3, see text for explanation. Parameters used: N=33, 
ly = 10-^^ 7) = 0.5 • 10-^ u' = 1 



evolution operator in S, leaves / invariant, that is T*[y(0)] = y{t) G /. The 
kinetic subspace K C S, defined by y{0) = (u„, 0) is obviously the usual fluid 
GOY model. Further subspaces are the Alfvenic subspaces A^ defined by 
y(0) = (w„,±w„), say Z+ ^ and Z~ = (or vice versa). Each initial con- 
dition in these subspaces is actually a fixed point of the system. We studied 
the properties of stability of K and A^ . Following p^ , let us define for each 
/ the orthogonal complement P, namely S — I (B P. Let us then decompose 
the solution as y(t) — {yint{t), yext{t)) where the subscripts refer to the / and 
P subspaces respectively. Finally we can define the energies Eint = \\yi7it\\'^ 
and Eext = ||yext|P- Note that the distance of a point y = {yint,yext) from 
the subspace / is d = min ||j/ — j/|| = ||yea;t||- Then Eext represents the square 

of the distance of the solution from the invariant subspace. At time t — 0, 
Eext = sEint (e <^ 1) represents the energy of the perturbation. Since the 
total energy is constant in the ideal case, two extreme situations can arise: 1) 
The external energy remains of the same order of its initial value, that is the 
solution is trapped near / which is then a stable subspace; 2) The external 
energy assumes values of the same order as the internal energy, that is the 
solution is repelled away from the subspace which is then unstable. Since the 
external and internal energies for the Alfvenic subspaces are nothing but the 
pseudocnergics E^ and E~ , which are ideal invariants, the Alfvenic subspaces 
are stable. As regards the kinetic subspace, Eint and Eext represent respec- 
tively the kinetic and magnetic energies. Looking at the numerical solutions 
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Fig. 3. Ideal case: kinetic energy {continuous line) and magnetic energy {dashed 
line) versus time for the 3D model; magnetic energy {dot-dashed line) versus time 
for the pseudo 3D model 

of the ideal model (fig. ||) we can see the difference in the stability properties 
between the pseudo 3D model and the 3D one. In the first case the external 
energy remains approximately constant, while in the second case the system 
fills up immediately all the available phase space. This striking difference is 
entirely due to the nonlinear term, and in fact must be ascribed to the differ- 
ences in sign of the third invariant. The effect of the unstable subspace, which 
pushes away the solutions, is what in rcf. p6[ is called "turbulent dynamo 
effect". 



5 Spectral properties in stationary forced state 

The main fundamental difference between hydrodynamic and MHD shell 
models lies in the fact that the behaviour of the former is not so sensitive to 
the type of forcing, at least as far as the main features are concerned. On the 
contrary in the magnetic case phase space is more complex because of the 
presence of invariant subspaces which can act as attractors of the dynamics 
of the system, hence the type of forcing becomes crucial in selecting the sta- 
tionary state reached by the system. The spectral properties of the 3D model 
have been investigated by Frick and Sokolov in [ [l6[ under different choices 
of the forcing terms. In their simulations they observe that the spectral in- 
dexes of kinetic and magnetic spectra depend on the level of cross helicity 
and magnetic helicity. In particular spectra with spectral index —5/3 appear 
if the cross helicity vanishes. Even in this case results may be deceptive. In 
fact, defining the reduced cross helicity hfj as the cross helicity divided by 
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the total energy, long runs flj show that, in case of constant forcing on the 
velocity variables, even from an initial value hjf — the system evolves in- 
evitably towards a state in which the reduced cross helicity reaches either the 
value +1 or -1, corresponding to a complete correlation or anti-correlation 
between velocity and magnetic field. In terms of attractors the system is at- 
tracted towards one of the Alfvenic subspaces where velocity and magnetic 
field are completely aligned or anti-aligned. Due to the particular form of the 
nonlinear interactions in MHD (0) , the nonlinear transfer of energy towards 
the small scales is stopped. In this case Kolmogorov-like spectra appear as 
a transient of the global evolution. This is shown in fi-g.^ where it is clearly 
seen a component (Z^) which is completely vanishing while the Z+ spectrum 
becomes steeper and steeper as energy is not removed from large scales. If 
an exponentially correlated in time gaussian random forcing on the velocity 
field is adopted the system shows a very interesting behaviour. It spends long 
periods (several large scale turnover times) around one of the Alfvenic attrac- 
tors, jumping from one to the other rather irregularly (fig. |5|). This behaviour 
assures the existence of a flux of energy to the small scales, modulates the 
level of nonlinear interactions and the consequent dissipation of energy at 
small scales, which is burstly distributed in time. What we want to stress is 
the fact that the Alfvenic attractors play a relevant role in the dynamics of 
the system. This fact should be taken into account especially when a station- 
ary state is investigated in order to determine the scaling exponents of the 
structure functions (see below). Two regimes, the Kolmogorov transient and 
the completely aligned regime, could be mixed during the average procedure, 
thus leading to unreliable values of the scaling exponents. 

6 Fluxes, Inertial Range and Intermittency 

The "four-fifth" relation {6v{l)^) ~ (— 4/5)e?, where e is the mean rate of 
energy dissipation and / the separation, derived by Kolmogorov in |23], can 
be generalized to MHD flows |23| , [g^ , |g^ . A corresponding relation exists in 
MHD shell models, which can be derived following the considerations in pfl . 
Assuming for simplicity v — rj, the scale-by-scale energy budget equation is: 

2—1,71 2 — 1,71 2— l,n 

where the quantities {Z^Z^Z^)^ are defined as 



r±Z±ZT)„ = llm{(<5 + S„,)}Z^Z^^,Z^^, + I^-^-i!lllz±_i^,T^„% 



+ {2-S- S„,)ZtZ^^,Zt+, + i^V^Z„T_iZ±Z±+i} (9) 



Assuming that i) forcing terms only act at large scales; (ii) the system tends 
to a statistically stationary state; (iii) in the infinite Reynolds numbers limit 



10 



Paolo Giuliani 



0- 



N 



60 4 




oooo 



OO 



OOoOOO 



o 



0- 



N 



bO 4" 




1 



2 3 4 5 
log,„ k„ 



12 3 4 5 6 7 
logio k„ 



Fig. 4. Left: logj^Q(|Z^|) versus logjQ fc„ at different times (a) Black circles: The 
average of \Z^\ are made over the first 30 large scale turnover times (b) White 
circles: The average is made after 300 large scale turnover times (c). The straight 
line has slope —1/3. Right: The same for \Z^\ 
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Fig. 5. Reduced cross helicity versus time in the case of an exponentially time- 
correlated random gaussian forcing on the velocity variables 
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{v — > 0) the mean energy dissipation tends to a finite positive limit e , we 
obtain 

{{Z+Z+Z-)„) = -e+k-' 

{{Z-Z-Z+)„) = -e-k-' 

These are the equations that define the inertial range of the system and that 
can be easily checked and confirmed by numerical simulations (fig. |g) . It is to 
be remarked that these are the appropriate combinations that are expected to 
scale exactly as fc~^. Let us finally remind that, as far as cascade properties 
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of shell models are concerned, the major drawback lies in the difficulty to 
reproduce cascades of quantities that are expected to flow inversely, such as 
energy in 2D hydrodynamic EG] or magnetic helicity in MHD H. 

A deep understanding of intermittency in turbulence is nowdays one of the 
most challenging tasks from a theoretical point of view (see [27) for review) . 
A lot of papers have been dedicated in the last years to investigate tempo- 
ral intermittency in shell models. Deviations from the Kolmogorov scaling 
?P = p/3 of the scaling exponents in the structure functions, (|Mn|)^ ~ k~^p , 
have been observed and described in the context of a multifractal approach 
[ p^ . A precise calculation of the scaling exponents may have difficulties re- 
lated to the presence of periodic oscillations superimposed to the power law. 
Another source of uncertainty is linked to the exact identification of the in- 
ertial range where the fit should be performed. These problems are at lenght 
discussed and investigated in p8] where a new shell model (called Sabra 
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model) has been introduced in the context of hydrodynamic turbulence. The 
Sabra model is a slight modification of the standard GOY model and allows 
to eliminate spurious oscillations in the spectra. The same problems are in 
principle encountered in magnetohydrodynamic models thus a generalization 
of the Sabra model to MHD is required (pol). An alternative approach to 
the determination of scaling exponents for the 3D MHD GOY model can be 
found in pQ] where concepts and techniques related to ESS and GESS [pl[ 
are used. 

We have determined the scaling exponents of the structure functions 
{\un\P) - A:-«p, {\bn\P) - fc-C, {\Z+\P) - fc-C, {\Z-\P) ^ fc-«p adopting a 
random forcing on the velocity variables (on shell n=l and ii—2) to assure 
the system does not "align" . The forcing terms were calculated solving a 
Langevin equation /„ = — (I/tq) /« + fJ., where tq is a correlation time chosen 
equal to the large scale turnover time and fi is a gaussian delta-correlated 
noise. The total number of shells is 23 and the values of viscosity and resis- 
tivity are v = 0.5 • 10~^, 77 — 0.5 • 10~^. In fig. ^ the first three structure 
functions are plotted for the magnetic field, together with the best fit lines. 
From a comparison with spectra obtained in the standard GOY model Ig^l j 
it should be remarked that the cross over region between the inertial range 
and the dissipative one is not so sharp as in the hydrodynamic case. 

We then decided to perform a least-square fit in the range, determined 
visually, between the shell numbers n = 3 and n — 12. The values of ^" and 
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^p are reported in Table |l| together with the values of ^p, extracted from |^ , 
for the hydrodynamic GOY model. The values of the scaling exponents of the 
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other structure functions are compatible, within errors coming from the fit 
procedure, with those of the velocity variables. It can be seen that the values 
found are compatible with those obtained for the standard hydrodynamic 
GOY model. 



Table 1. Scaling exponents ^p (GOY model), ^p, ^^ 



p 


Cp (GOY) 




Sp 


1 




0.37 ±0.01 


0.36 ±0.01 


2 


0.72 


0.71 ±0.01 


0.70 ±0.01 


3 


1.04 


1.02 ±0.02 


1.02 ±0.02 


4 


1.34 


1.31 ±0.03 


1.32 ±0.03 


5 


1.61 


1.57 ±0.04 


1.59 ±0.04 


6 


1.86 


1.82 ±0.05 


1.84 ±0.05 


7 


2.11 


2.05 ±0.06 


2.07 ±0.06 


8 


2.32 


2.28 ±0.08 


2.30 ±0.08 



7 Conclusions 



In this paper we have reported about the main properties concerning dynam- 
ical behaviour and intermittency of a shell model for MHD turbulence. The 
properties of the model reveal a complex structure of phase space in which 
invariant subspaces are present. The stability properties of the kinetic sub- 
space are related to a dynamo action in the system while Alfvenic subspaces 
act as strong attractors which cause the system to evolve towards a state in 
which no energy cascade is present. A careful choice of forcing terms seems 
to be crucial in determining the stationary state reached by the system. 

We want finally mention that shell models, as good candidates to repro- 
duce the main features of MHD turbulence, can be used to check conjectures 
and ideas in astrophysical applications where very high Reynolds numbers are 
often present. We briefly remind two examples of applications. In p3] MHD 
shell models have been used to simulate magnetohydrodynamics in the early 
universe to investigate the effects of plasma viscosity on primordial magnetic 
fields. As second example, scaling laws found in the probability distribution 
functions of quantities connected with solar flares (eruption events in the so- 
lar corona) are at present matter of investigation by means of shell models. 
Results on this subject can be found in pJ. 
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